Similarity of fluctuations in correlated systems: The case of seismicity* 
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We report a similarity of fluctuations in equilibrium critical phenomena and non-equilibrium 
systems, which is based on the concept of natural time. The world-wide seismicity as well as 
that of San Andreas fault system and Japan are analyzed. An order parameter is chosen and its 
fluctuations relative to the standard deviation of the distribution are studied. We find that the scaled 
distributions fall on the same curve, which interestingly exhibits, over four orders of magnitude, 
features similar to those in several equilibrium critical phenomena ( e.g., 2D Ising model) as well as 
in non-equilibrium systems (e.g., 3D turbulent flow). 
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I. INTRODUCTION 

Recently, a great interest has been focused on the fluc- 
tuations of correlated systems in general and of critical 
systems in particular0,13, ^, 4,5, .6, 8, _9J. Bramwell, 
Holdsworth and Pinton (BHP)[y, in an experiment of a 
closed turbulent flow, found that the (normalized) prob- 
ability distribution function (PDF) of the power fluctu- 
ations has the same functional form as that of the mag- 
netization (M) of the finite-size 2D (two-dimensional) 
XY equilibrium model in the critical region below the 
Kosterlitz-Thouless transition temperature (Magnetic or- 
dering is then described by the order parameter M). The 
normalized PDF, denoted by P{rn), is defined by intro- 
ducing the reduced magnetization "li m = (M — {AI})/a, 
where (M) denotes the mean and a the standard devia- 
tion. For both systems, BHP found that while the high 
end (to > 0) of the distribution hasQ a Gaussian shape 
the asymptote of which was later clarified^ to have a 
double exponential form, a distinctive exponential tail 
appears towards the low end (to < 0) of the distribution. 
The latter tail, which will be hereafter simply called, for 
the sake of convenience, "exponential tail" , provides the 
main region of interest |l|, since it shows that the prob- 
ability for a rare fluctuation, e.g., of greater than six 
standard deviations from the mean, is almost five orders 
of magnitude higher than in the Gaussian case. Subse- 
quent independent simulations 0, IE S IE showed that 
a variety of highly correlated (non equilibrium as well as 
equilibrium) systems, under certain conditions, exhibit 
approximately the "exponential tail" . 

Earthquakes do exhibit c omp lex correlations in space, 
time and magnitude, e.g. flOl lTll IT^. It has been re- 
peatedly proposed (see Ref.[i3] and references therein) 
that the occurrence of earthquakes (cf. mainshocks) 
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can be considered as a critical point (second-order phase 
change), but alternative models based on first-order 
phase transitions have been also forwarded which are 
probably more applicable, see Ref.|0| and references 
therein. (Such a diversity also exists for the brittle rup- 
ture which is a phenomenon closely related to earth- 
quakes. Buchel and SethnaJ5|| have associated brittle 
rupture with a first-order transition and a similar view 
has been also expressed in Re fs.jT^ IT^ . On the other 
hand, Gluzman and Sornette[lE| later suggested that 
it is analogous to a critical point phenomenon.) Both 
approaches lead to scaling laws or power-law distribu- 
tions for the dynamical variables (second-order transition 
demonstrate scaling near a critical point, whereas first- 
order transitions demonstrate scaling when the range of 
interactions is large (mean-field condition), as is the case 
with elastic interactionsHj). However, the question on 
whether earthquakes exhibit an "exponential tail" , has 
not yet been clarified. This might be due to the major 
difficulty of choosing an order parameter in the case of 
earthquakes (EQs). Following the wording of Ref . [l^ . we 
note that in general such a choice is an art, since usu- 
ally it's a new phase which we do not understand yet, 
and guessing the order parameter is a piece of figuring 
out what's going on. The scope of the present paper 
is twofold: to propose an order parameter for the case 
of EQs and then examine whether an "exponential tail" 
appears. We find that our scope is achieved only if we 
analyze the series of earthquakes in the natural time- 
domain HE HI HI HI H El El 113 . 

In order to serve the aforementioned scope, the present 
paper is organized as follows: In Section II, we explain 
how the power spectrum of the seismicity in natural time 
can be obtained. An order parameter for EQs is proposed 
in Section III. In the light of this proposal, and without 
using any adjustable parameter, we show in Section IV 
that the normalized distribution of the long term seismic- 
ity for different seismic areas fall on a universal curve. It 
consists of two segments the one of which exhibits the 
"exponential tail" . Interestingly, a further investigation 
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of the latter segment in Section V, reveals that it is sim- 
ilar to that observed in several equilibrium critical phe- 
nomena (e.g., 2D Ising, 3D Ising) and in non-equilibrium 
systems (e.g., 3D turbulent flow). A brief discussion fol- 
lows in Section VI, while Section VII summarizes our 
main conclusions. Two Appendices provide clarifications 
on some points discussed in the main text. 

II. THE SEISMICITY IN NATURAL TIME 

In a time series consisting of N events, the natural time 
Xk = k/N serves as an index pol l2l| for the occurrence 
of the fc-th event. It is, therefore, smaller than, or equal 
to, unity. For the analysis of scismicity, the evolution of 
the pair {xk,Ek) is considered 20„. .22> .30. .31.] . where Ek 
denotes the seismic energy released during the k-th event 
see FigQ] (cf. This energy -which is itself proportional to 
the seismic moment Mq and hence we can use in the verti- 
cal axis of Fig^b) either Ek or (Afo)/c is related to the 
magnitude M through E (x where c is a constant 

around 1.5). The following continuous function F(uj) 

was introducedlMIIlIll: F{uj) = Ef=i ^fc exp (iw|:) 
where uj = 27r(/), and (f) stands for the natural frequency. 
We normalize F{uj) by dividing it by F(0), 

Z^n=l fe=l ^ 

where pk = En- A kind of normalized power 

spectrum n(a;) can now be defined: n(w) = |$(ti;)|^. 

For a Seismic Electric Signals (SES) activity, which is 
a sequence of low frequency (< 1 Hz) electric pulses emit- 
ted when the stress in the focal area approacheS i28i. i32l 
a critical value, we have shown (for details see Ref . 1201 . 
see also [IJ) that the following relation holds[2(ill2lll22l| 



III. THE ORDER PARAMETER PROPOSED 

We now proceed to choose the order parameter, as- 
suming that a mainshock may be considered as the new 
phase. We take advantage of the experimental factp3| 
that several hours to a few months before a mainshock 
an SES activity is recorded, and focus our attention on 
the evolution of the seismicity (in the candidate area) 
during the period from the SES detection until the main- 
shock. If we set the natural time for the seismicity zero 
at the initiation of the concerned SES activity, we form 
time series of seismic events in natural time (see Fig^ 
for various time windows as the number N of consecutive 
(small) EQs increases. When computing Il{(f>) (as well as 
Ki, see below) for each of the time windows, we find that, 
in the range < < 0.5, it approaches, as N increases 
from 6 to some value less than (or equal to) 40, to that 
given by Eq.® (or the Ki-value becomes equal to 0.07, 
see Eq.Q). The coincidence occurs only a few hours to 
a few days before the mainshock. (In simple words, be- 
fore a mainshock a sequence of earthquakes occurs, which 
obeys Eq.(|2J) and this process will be called single cor- 
related process.) When the mainshock occurs (the new 
phase), n(0) abruptly increases to approximately unity 
(for details see Ref. . 29j ) and ki becomes almost zero. 
This can be visualized in the example depicted in Fig|21 
where we plot the Ki-value versus the number of EQs af- 
ter the SES detection on April 18, 1995 (see Refs.^1, 28]) 
until the occurrence of the M=6.6 mainshock on May 
13, 1995 at 'i0.2°N,21.7°E. This figure shows that the 
Ki-value becomes ki w 0.07 after the 11th EQ (see also 
Ref.;295), while upon the mainshock the Ki-value abruptly 
decreases to ki « 9 x lO"'^. Such a behavior has been 
verified[23, H^l for several major EQs and points to the 
conclusion that Tl{(f>) for small 0, or ki, could be consid- 
ered as an order parameter. 
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We focus on the properties of n(cj) or 11(0) for natural 
frequencies cj) less than 0.5, since in this range of (p, n(ijj) 
or n((/>) reduces I20I I21L [22. 1301 to a characteristic func- 
tion for the probability distribution pk in the context of 
probability theory. According to the probability theory, 
the moments of a distribution and hence the distribution 
itself can be approximately determined once the behavior 
of the characteristic function of the distribution is known 
around zero. For uj 0, Eq.© leads to[20ll2lll3(il|: 



U{lo) « 1 - 0.07^2 



(3) 



which reflects (see Appendix that the variance of x is 
given by: 



(4) 



In Section IV, we will investigate whether Eq.(|2Jl holds 
for EQs. 



IV. UNIVERSAL CURVE FOR SEISMICITY 

The properties of the power spectrum for the long term 
seismicities in natural time can be studied by means of 
the following procedure: First, calculation of 11(0) was 
made for an event taking time windows from 6 to 40 con- 
secutive events (for the reasons explained in Section III; 
the choice of the precise value of the upper limit, up to 
100 or so, is not found decisive |23, 122, ISJ)- And second, 
this process was performed for all the events by scanning 
the whole catalogue. The following data from two dif- 
ferent areas, i.e., San Andreas fault system and Japan, 
have been analyzed: First, the EQs that occurred dur- 
ing the period 1973-2003 within the area N^^W^^i us- 
ing the Southern California Earthquake catalogue (here- 
after called SCEC). Second, the EQs within 7V||^^i25 
for the period 1967-2003 using the Japan Meteorologi- 
cal Agency catalogue (hereafter simply called "Japan"). 
The thresholds M>2.0 and M>3.5 have been considered 
for SCEC and Japan, respectively, for the sake of data 
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completeness 29]. By plotting for a given value of 4> the 
observed probability P[n((/))] versus !!((/)) (two such ex- 
amples are given in FigsOfa) andOl^b) for cf) = 0.05 and 
(j) = 0.4, respectively), we find that a local maximum 
occurs at a value of n((/)) hereafter called Ii-p{4>) (see 
also Appendix EJ- This lies very close (see FigOJc)) to 
the value I{th{4>) obtained theoretically, i.e., estimated 
from Eq. Q . The validity of Eq. Q for various 0- values, 
in the range < cf) < 0.5, can be now visualized in 
FigEId), where we see that np(0)-values versus <j) for 
both SCEC and Japan do not differ by more than 1% 
from the lit ^ (</))- values (cf. this difference is more or less 
comparable to the estimation error of np((/)), for details 
see Appendix Ib]|. 

We now plot, in Fig^ the quantity (tP{X) versus 
{X - {X))/a where X stands for n((/)) and (n(0)) and 
a refer to the mean value and the standard deviation 
of n(0) (recall that the calculations should be done for 
small 0- values, e.g., (/)=0.05, since we assume here — » 0, 
for the reasons explained in Section II) . One could alter- 
natively plot aniPini) versus ((ki) — ki)/ctki, where (ki) 
and CT„j now refer to the mean value and the standard 
deviation of ki. The results in Fig. 21 for both areas, 
fall on the same curve. This log-linear plot clearly con- 
sists of two segments: The segment to the left shows a 
decrease of P{X) almost by five orders of magnitude, 
while the upper right segment has an almost constant 
P{X) (Obviously, the latter segment deviates from the 
general behavior of the BHP distribution -as it was sum- 
marized in Section I- but from thereon we put emphasis 
on the left segment since our main interest here concerns 
the "exponential tail"). The feature of this plot is strik- 
ingly reminiscent of the one obtained by Bak et al.jlOj 
(see their Fig. 4) on different grounds, using EQs in 
California only. More precisely, they measured Ps\i{T), 
the distribution of waiting times T, between EQs occur- 
ring within range / whose magnitudes are greater than 
M = \ogS. They then plotted T°'Ps,i{T) versus TS-H'^ 
and found that, for a suitable choice of the exponent a 
(i.e., a = 1), the Gutenberg-Richtcr law exponent h (i.e., 
h=l) and the spatial dimension d (i.e., fractal dimension 
d = 1.2) all the data collapse onto a single curve which 
is similar to that of Fig0J Recall, however, that Fig0] 
was obtained here without considering at all the wait- 
ing time distribution and without the suitable choice of 
any parameter. After a further inspection of Fig^ the 
following points have been clarified: 

First, the rapidly decaying part (i.e., the left segment), 
which is consistent with an almost exponential decay- 
ing function over almost four orders of magnitude, re- 
mains practically unchanged, upon randomizing the data 
("shuffling" 123). (Some changes do occur in the right 
part, associated with aftershocks, see also below.) This 
can be seen in the inset of Fig. ^ where for the sake 
of clarity only the results from the data of Japan (the 
original as well as the "shuffled" ones) are depicted. 

Second, the feature of the plot of Fig. 0] is not al- 
tered upon changing either the seismic region or the time- 



period (provided that the latter does not include after- 
shocks on/j/, see below). As an example, FiglSJb) shows 
that three different regions A, B, C in Japan (depicted 
in Fig[5Ja)), as well as the whole Japan, result in almost 
identical plots. 

Third, the "upturn branch" in the upper right part of 
Fig. 21 arises from the presence of aftershocks. It dis- 
appears (see the crosses in Fig. when, in Japan, for 
example, we delete the EQs with M<5.7 (and hence dras- 
tically reduce the number of aftershocks), but it does not, 
when deleting EQs with smaller threshold, i.e., M<4.0; 
the latter can be also visualized in the SCEC example 
of Fig. IHl where we give the results for M>4.0 (cf. this 
threshold still allows the presence of a reasonable number 
of aftershocks). 

Fourth, if we consider the relevant results for the 
worldwide seismicity (WWS) by taking a large magni- 
tude threshold, i.e., M>5.7 (so that for the data to be 
complete pgjl'). we find (see Fig|Hlthat will be further dis- 
cussed below) that they fall onto the same curve with the 
results of both Japan and SCEC. 



V. DOES A UNIVERSAL BEHAVIOUR EXIST 
FOR DIVERSE SYSTEMS? 

We now compare in FigEl the aforementioned re- 
sults of seismicity with those obtained in some equi- 
librium critical systems (e.g., see Ref.Q)- We first re- 
call that the PDF in the critical regime depends on 
K — 1/T and the length L through a scaling variable 
s = L^/^iK - Kc)/Kc, where = l/T^ and de- 
notes the critical temperature (the quantity provides 
the ratio of the lattice size and the correlation length 
at K). In Fig. El we include numerical results of the 
2D Ising model for s = 8.72(L = 128, K = 0.4707) and 
s = 17.44(L = 256, K = 0.4707). Here, X stands for M. 
These s-values were intentionally selected, because 8] for 
s > 8.72 for the 2D Ising model, the P{m, s)'s of a num- 
ber of critical models (i.e., 2D XY, 2D Ising, 3D Ising, 2D 
three-state Potts) share the same form (up to a constant 
factor of s), which interestingly exhibits an exponential- 
like left-tail (to < 0). An inspection of Fig|Hl shows that 
our 2D Ising results almost coincide (cf. this can be safely 
checked only for the left segment, i.e., to < 0) with those 
of seismicity, i.e., Japan, SCEC and WWS (cf. Some dis- 
parity which appears in the upper right part of SCEC 
only, might be attributed to the selection of the magni- 
tude threshold for seismicity, recall the third point men- 
tioned in Section IV). This coincidence (which seems to 
be better for s — 17.44) reveals that the seismicity, ir- 
respective of the seismic area we consider, exhibits -over 
four orders of magnitude- fluctuations of the order pa- 
rameter similar to those in several critical systems as well 
as in 3D turbulent flow. 
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VI. DISCUSSION 

It is of interest to see how the scaled distributions look 
like in the frame of the present analysis, if one generates 
surrogate data either by means of a simple Poisson model 
or by the Gutenberg- Richter law and compare the results 
to those deduced from actual seismicity data. 

In Fig[7| we present the linear-linear plot (Fig. EJa)) 
as well as the log-linear plot (Fig. db)) of aP{X) ver- 
sus {X — {X))/a where X stands for 11(0) for « 0, 
for surrogate data of EQs for which their {Mo)k obey a 
simple Poisson rule for various mean values A lying be- 
tween 5 and 200. In the same figure we insert the results 
for Japan (M > 3.5) and SCEC (M > 2.0) aheady dis- 
cussed in Fig^ Although we find that upon decreasing 
A the surrogate data move closer to the real data, how- 
ever a satisfactory agreement between them cannot be 
supported. 

In Fig. IHlwe repeat the procedure followed in Fig. [7| 
but now the surrogate data are produced on the basis 
of the Gutenberg-Richter law, i.e., that the (cumulative) 
number of EQs with magnitude greater than M (occur- 
ring in a specified region and time) is given by 



N{> M) « 10 



-6M 



(5) 



It is currently considered^] that b is generally a con- 
stant varying only slightly from region to region being 
approximately in the range 0.8 < b < 1.2. For Japan 
and SCEC we find on the basis of Eq.@ b w 1.05. Note 
that in FiglHl surrogate data are intentionally produced 
for a variety of b values in the range b = 0.5 to 2.0. 
An inspection of this figure leads to the following con- 
clusions: First, the curves of the surrogate data marked 
with b = 0.5 tob = 0.9 significantly differ from that of the 
real data. Second, for &- values larger than 1 and smaller 
than 1.4, the curves of the surrogate data have a general 
feature more or less similar to the curve of the real data. 
However, none of these &- values in the surrogate data can 
lead to a curve coinciding to the one obtained from the 
real data. 

In other words, the scaled distribution, deduced within 
the frame of the present analysis, reveals for the real 
data an extra complexity when compared to the surro- 
gate data even if the latter are produced with 6-values 
comparable to the experimental ones. 



VII. CONCLUSIONS 

The main conclusions could be summarized as follows: 

(1) The analysis of the seismicity in the natural time- 
domain reveals that Ii{4>) (for small 0) or ki, may be 
considered as an order parameter. 

(2) If we study the order parameter fiuctuations relative 
to the standard deviation of its distribution, the following 
two facts emerge (without making use of any adjustable 
parameter) : 



First, the scaled distributions of different seismic areas 
(as well as that of the world wide seismicity) fall on the 
same curve (universal) . 

Second, this curve exhibits an "exponential tail" form 
similar to that observed in certain non-equilibrium sys- 
tems (e.g., 3D turbulent flow) as well as in several(e.g., 
2D Ising, 3D Ising, 2D XY) equilibrium critical phenom- 
ena. 



APPENDIX A: DERIVATION OF EQUATION (H) 

The Taylor expansion, around u;=0, of the relation 
n(w) = \<^{ui)f using Eq(l) reveals that [13 

Il{uj) = 1 - KlUj"^ + K2LU^ + K3LJ^ + + ■ ■ ■ (Al) 

where 

1 d^n{uj) 



2 dw^ 



(A2) 



j=0 



We now consider 



^ -$*(w) — -\--+<^{uj) — ^^^-f2- ^ ' ^ ' 
dui 



div^ 



diJ^ 



duj dio 
(A3) 

and taking into account that $(0;) = exp(iti;xfc), 
with $(0) = 1, we find: 




= (x^)-(x)^ 

where (x") = J2kPkXk- 

Expanding Eq.|(2Jl, around lo — 0, we have: 

n(w) = 1 - 0.07^2 -I- . . . 

and hence: 

= ix') - {X? = 0.07 
which is just Eq.lQ. 



(A4) 

(A5) 
(A6) 



APPENDIX B: THE PROCEDURE TO 
DETERMINE THE MAXIMUM IN P[n(<?!))] 
VERSUS n(0) 

The calculation of 11(0) was made, as mentioned in 
Section IV, for an event taking time windows from 6 to 
40 consecutive events and this process was performed for 
all the events by scanning the whole catalogue. This 
procedure resulted in the calculation, for each 0-value 
and each catalogue, of more than 10® n(0)-values, whose 
probability density function (PDF) was determined by 
using the computer code histogram of Ref . |34j with a 
number of bins proportional to iV^^'^, where N is the 
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number of n(0)-values (cf. this point, i.e., that the num- 
ber of bins should be proportional to N^^^, is discussed 
in Ref. 1^3)- This method resulted in the PDFs shown 
in Figl^l as well as in those depicted in Figs|3Ja)|2Ib), 
01 Elb) , and |S1 Due to the intrinsic fluctuations of the 
values of the PDF (because N is still finite), a direct 
determination of the value np(0) where the PDF maxi- 
mizes, just by simply taking the maximum value of the 
calculated PDF, may lead to erroneous values of 11^(0). 
One should consider the general trend of the PDF as a 
whole, which can definitely show a more accurate and 
stable value of Ilp{4>). Thus, the procedure we applied 
for the determination of Ilp{(j)) was as follows: For each 
(/)- value, a region [a, b] around the maximum was se- 
lected (examples are shown in FigsEIa) and (b)) and 
then a cubic polynomial, p{x) — a + bx + cx^ + dx^, 
was used to fit the PDF values in this region. (Close 



enough to the maximum, a parabolic fit could be also 

good since f{x) = frnax - \fmax\{^ " Xmaxfl'^, but 

in view of the PDF asymmetry the cubic polynomial 
used, provides a better approximation in the whole re- 
gion [a, b].) The value of T\p[(j)) was determined through 
the direct maximization of this cubic polynomial, i.e, 

FiglSJc) have been obtained by means of such a proce- 
dure. Finally, we note that, due to the fitting procedure 
involved and the relative arbitrariness in the definition of 
[a,b], the estimation error of Wp{(j)) is more or less com- 
parable to its percentage deviation from Tith{4')^ depicted 
in FiglS^d). Thus, we can state that np(0) and Ilth{4i) 
are experimentally indistinguishable, which strengthens 
the statement that Eq.(2) -which has been used for esti- 
mating I{th{<j>)- holds for EQs. 
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FIG. 1; How a series of seismic events in conventional time 
(a) is read in the natural time (b). This example refers to 
the first If small earthquakes (cf. the month/date is marked 
on the horizontal axis in (a)) that occurred after the SES ac- 
tivity recorded on April f8,f995 and preceded the mainshock 
(M=6.6) of May f3,f995. 
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FIG. 2: How the variance ki evolves event by event during 
the following period: from the detection|28| of the SES ac- 
tivity on April f8, f995 until the occurrence of the M=6.6 
mainshock(labelled f8) on May f3, f995. All the EQs used 
in the calculation are tabulated in Ref. 29]. (cf. the first fl 
EQs -out of fS- are those depicted in FigQ^a)). 
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FIG. 3: (Color online) Validity of Eq.@ for SCEC and 
Japan. We first determine (see also Appendix B) for 
SCEC(circles) and Japan(crosses), for each value, the value 
^p{4>) at which P[n((/>)] maximizes. Two such examples are 
shown in (a) and (b) for (j> = 0.05 and (j> = 0.4, respec- 
tively. In (c), we plot the resulting 11^(0) values versus (p 
for SCEC(circles) and Japan(crosses); the solid line corre- 
sponds to Eq.(|5J. The percentage difference between Ilp{4') 
and nth(<^) (obtained from Eq.(|5J) is shown in (d). 
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FIG. 4: (Color online) Universality of the probability density 
function of Tl{4>) for EQs in the natural time-domain. The 
log- linear plot of aP{X) versus {X — {X))/a, where X stands 
for n((j!>) for 0^0. Crosses and circles correspond to Japan 
(M>3.5) and SCEC (M>2.0), respectively. The inset depicts 
the corresponding results for the "shuffled" data (black curve) 
and the original data (red crosses) in Japan. The same graph 
is obtained for three different regions in Japan (see Fig|SJ. 




FIG. 5: (Color online) The same as Fig. 0] but for the regions 
A (red), B (green) and C (blue) in Japan (b). A map of these 
regions is shown in (a). 
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FIG. 6: (Color online) The common feature of fluctuations 
in different correlated systems. The log-linear plot of aP{X) 
versus {X — {X))/a for WWS (triangles), Japan(crosses) and 
SCEC(circles). The magnitude threshold M>5.7 for WWS 
and Japan (while M>4.0 for SCEC) was used see the text. 
Furthermore, the dotted curve shows the results obtained for 
the 2D XY model fwith[3^ inverse Kosterlitz-Thouless tran- 
sition temperature Kkt ~ 1.2) {X = ^M| + M^), K = 2.0 
for L = 10(A'^ = 100) which has been shownU to describe the 
experimental results for 3D turbulent flow. The results of the 
2D Ising model K = 0.4707 (while Kc « 0.4407), either for 
L = 128(dashed) or L = 256(solid line), are also plotted. 
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FIG. 7: (Color online) The linear-linear (a) and log-linear 
(b) plots of aP{X) versus {X — {X))/a, where X stands for 
n(<j!>) for « 0. The results for Japan (M> 3.5) and SCEC 
(M> 2.0) are plotted along with those deduced from a series of 
independent and identicaly distributed (Mo)*: sampled from 
a Poisson distribution with mean value 5, 10, 20, 30, 40, 50, 
60, 70, 80, 90, 100 and 200: from the bottom to the top in the 
maxima appearing in (a) , respectively, and from the upper to 
the lower left branch in (b), respectively. 
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FIG. 8: (Color online) The linear-linear (a) and log-linear 
(b) plots of aP{X) versus {X — {X))/a, where X stands for 
n(<j!>) for « 0. The resuhs for Japan (M > 3.5) and SCEC 
(M > 2.0) are plotted along with those deduced from shuffled 
artificially generated EQ data obeying the Gutenberg-Richter 
law for various values of the exponent b =0.5, 0.7, 0.9, 1.0, 
1.1, 1.2, 1.3, 1.5 and 2.0 from the lower to the upper curve at 
the value {X — {X))/a ~ —0.5 in (a), and from the upper to 
the lower curve at the value [X — {X))/a ~ —4 in (b). 
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FIG. 9: (Color online) Two examples of the procedure used 
in the determination of Hp ((/)). (a) For SCEC at <^ = 0.03, (b) 
For Japan at = 0.4. The (red) curves in each case, show 
the cubic polynomial fit which was used in the range [a, b] 
around the maximum. The (red) arrow indicates the position 
of np(0). 



